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Abstract: We focus on the issue of proper definition of entanglement entropy in lattice 
gauge theories, and examine a naive definition where gauge invariant states are viewed as 
elements of an extended Hilbert space which contains gauge non-invariant states as well. 
Working in the extended Hilbert space, we can define entanglement entropy associated with 
an arbitrary subset of links, not only for abelian but also for non-abelian theories. We then 
derive the associated replica formula. We also discuss the issue of gauge invariance of the 
entanglement entropy. In the Z jv gauge theories in arbitrary space dimensions, we show 
that all the standard properties of the entanglement entropy, e.g. the strong subadditivity, 
hold in our definition. We study the entanglement entropy for special states, including the 
topological states for the Zn gauge theories in arbitrary dimensions. We discuss relations 
of our definition to other proposals. 
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1 Introduction 

Entanglement entropy plays important roles in various fields of quantum physics including 
string theory[l-8], condensed matter physics [9-14], and the physics of the black hole [15— 
19]. It is believed that entanglement entropy characterizes various aspects of quantum 
states in a simple and unified manner. 

In the context of lattice gauge theories, entanglement entropy is expected to be a useful 
tool for studying confinement / deconfinement transitions (or crossover) [20-22]. It has 
been pointed out, however, that there is a subtle problem in the definition of entanglement 
entropy in gauge theories[23, 24, 26, 27]. When we calculate the entanglement entropy of 
a region V, we first express the Hilbert space of the total system as a tensor product of the 
Hilbert spaces of V and that of V, the complement of V. Thus we trace out the degrees 
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of freedom of V and obtain the reduced density matrix of V. For gauge theories, however, 
the physical gauge invariant Hilbert space can not be factorized into a tensor product of 
the gauge invariant subspaces of V and that of V due to the local gauge invariance at 
the boundary dV between V and V. This reflects the fact that the fundamental physical 
degrees of freedom contain Wilson loops, which are nonlocal operators. Due to the absence 
of the factorization into a tensor product, it is not straightforward to define the reduced 
density matrix of some region and to calculate the entanglement entropy. We need to 
specify the prescription to obtain the reduced density matrix of the region. 

In this paper, we propose a definition of the entanglement entropy in lattice gauge 
theories. We extend the gauge invariant Hilbert space to a larger Hilbert space in order 
to admit the factorization into a tensor product of the gauge invariant subspaces of the 
region V and the region V in this larger Hilbert space. The natural candidate of this larger 
Hilbert space is the whole (gauge non-invariant) Hilbert space of the link variables. We 
then obtain the reduced density matrix of the region V by tracing out the link variables 
of the region V. We define the entanglement entropy as the von Neumann entropy of the 
above reduced density matrix. We can define the entanglement entropy for an arbitrary 
subset of links. This definition is applicable not only for abelian theories but also for non- 
abelian ones. We then derive the replica formula to calculate the entanglement entropy in 
our definition. 

In the gauge theories in arbitrary space dimensions, we express the whole Hilbert 
space by useful basis states, which are eigenstates of the gauge transformations [24]. We 
argue that all the standard properties of entanglement entropy, e.g. the strong subaddi¬ 
tivity, hold in our definition. We study the one for some special states. In particular, we 
calculate rigorously the one for the topological states in arbitrary space dimensions. We 
discuss relations of our definition to other proposals. We also demonstrate that the en¬ 
tanglement entropy depends on the choice of the gauge fixing for some simple cases. This 
indicates that one should not fix the gauge, at least on the boundary points between two 
regions, to calculate the entanglement entropy in gauge theories. 

The present paper is organized as follows. In section 2, we give precise definitions 
of the geometry of our lattice and define the entanglement entropy. We discuss the gauge 
invariance of the reduced density matrix. We also derive the replica formula here. In section 
3, we consider the Zn gauge theories. We express the whole Hilbert space by eigenstates of 
the gauge transformations, and derive an explicit expression of the entanglement entropy. 
We then argue that all the standard properties of entanglement entropy, e.g. the strong 
subadditivity, hold in our definition. We study the one for some special states. In particular, 
we calculate the one for the topological states in arbitrary space dimensions. We discuss 
relations of our definition to other proposals. In section 4, we summarize our investigations. 
Some properties of the Zj\j gauge theories used in the main text are given in appendix A, 
while gauge invariant states in non-abelian gauge theories are briefly discussed in B. 
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2 Naive definition of entanglement entropy in lattice gauge theories 
2.1 Definition and some properties 

Geometry We can treat quite general geometries and boundary conditions. 

Our lattice is ( S , £), where S denotes the set of sites x, y ,... € S, and C C S x S the 
set of links. We understand that £ = (x, y ) and £ = (y, x ) are different ways of expressing a 
single link in C. This in particular means that £ € V implies £ G V for any subset V C C. 
Here we do not assume a particular structure of our lattice such as regularity, so that a 
random lattice could be treated. Note that this setup can treat both periodic and free 
boundary conditions for the whole lattice. 

We define the boundary of a subset V C C as 

dV := { x € S | (x, y) € V and (x, z) G V for some y, z £ S }, (2.1) 

which is the set of sites in both V and its complement V = C\V. Note also that dV = dV. 

Naive definition of entanglement entropy We consider the global density matrix p 
for gauge theories, whose elements are denoted by 

(U\p\U') = p(U ; U') = p(Uv, Uy ; V v , U' v ) (2.2) 

where U represents a gauge configuration (a set of all link variables), U = {Ui \ £ e £}, while 
Uy, Uy are gauge configurations on V and V, Uy = {Up \ £ G V} and Uy = {Ug \ £ G V}, 
respectively. 

We propose to define a reduced density matrix as 

(Uy\pv\U!y) = py(Uy-U'y) = f VUy p(Uy ,Uy,U'y ,Uy) , (2.3) 

where T>Uy denotes a product of the group invariant integrals or sums. For the compact 
group, we have VUy = dUt, where dUi is the Haar measure for the link variable Ug. 

The above definition of the reduced density matrix is a simple generalization of the 
reduced density matrix in spin systems, where the whole Hilbert space is a direct product 
of those of region V and region V, Ji = Jiy <8> 'Hy- In the case of gauge theories, on the 
other hand, due to the local gauge invariance, the gauge invariant full Hilbert space can not 
be factorized into a product of gauge invariant subspaces, W^ v / Wy v <8> 'H“ v . Therefore 
the above reduced density matrix py can not be obtained from a single partial trace of 
p over the gauge invariant subspace "H^ v . Without gauge invariance, however, the whole 
Hilbert space can be factorized as He = Hy <S> Hy, so that our definition of py above 
can be understood as the partial trace of p over the gauge non-invariant subspace Hy on 
Uy. In the next section, we explicitly construct the reduced density matrix for the Z]y 
gauge theories in an arbitrary dimensions, and explicitly construct an extension of py to 
He = Hy <8> Hy. 

From the reduced density matrix, the entanglement entropy can thus be defined as 

S(V) = -tr \p v log p v ], (2.4) 
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where the trace is taken over Hy. The definitions (2.3) and (2.4) are so simple that they 
can be used not only for discrete abelian theories but also for continuous non-abelian gauge 
theories without practical difficulties. 

In the next section, we will see that this trace is reduced to a sum of traces in the gauge 
invariant subspace 'Hy v C Hy and discuss that a basic properties such as the symmetric 
property and the strong subadditivity are satisfied for the Z^ gauge theories. 

Gauge invariance In gauge theories, the global density matrix is gauge invariant as 

p(U° ] U' h )=p(U-,U'), (2.5) 

where the gauge transformation of the link variable Un is given by Uf = g x U^gl with 
t = (x,y). 

On the other hand, the reduced density matrix py in (2.3) does not have such a gauge 
invariance. Indeed, 

py{U^-U' v hv ) = J VUy P (Up,Uy-U^ hv ,Uy) = f VU^p{U^ , t/f; u!y hv , Uf) 

= j VUy p(Uy , Uy, {U!y) h ^v\Uy) = PV {U V] (E^)W), (2.6) 

where the invariance of the measure VU 9 ) = VUy and the gauge invariance of the full 

density matrix p such as 

p{U 9 v V ,U 9 y-U!y 9V ,U 9 y) = p(U V ,Uy-U!y,Uy) (2.7) 

are used. Therefore py is invariant under diagonal gauge transformations (gy = hy) only. 

This suggests that the reduced density matrix and thus its entanglement entropy may 
depend on the choice of the gauge if the gauge fixing is employed in the calculation. Indeed, 
we will show in the next section that values of the entanglement entropy are different for 
different gauge fixing conditions in some simple cases for the Zjy gauge theories. Because 
of this problem, it is important and sensible to calculate the entanglement entropy in the 
gauge invariant way without gauge fixing. 

2.2 Replica formula 

We briefly consider the replica formula for the entanglement entropy of lattice gauge the¬ 
ories based on our definition. 

Transfer matrix and path integral In lattice gauge theories, the evolution in a discrete 
time is given by the transfer matrix T (for example, see Refs. [28, 29]), which is given by 

T(U, U') := (U\T\U') = exp[^5 d (t/)]exp[5 0 (t/,t/ / )]exp[^ (i (t/ / )] (2.8) 
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where 


S#) = iE E Tr [V*,„Vr +i ,,uUl + ^Ui- 

” x£S fi=l,v=\ 

v¥=v 

d 


d 


2 r 


E E ^M 2 - 9 ) 


xSiS /i=l,i/=l 


So(U,U') = iEE Tr + 


25 


irEtS /i=l 


( 2 . 10 ) 


for the plaquette action on a ci-dimensional hyper-cubic lattice with the coupling constant 
g. We here define P fiv (x) as a trace of the plaqiutte on gu plane at x, and use the notation 
that U Xifl := Ufi with t = (x,x + p) and fi is an unit vector in the fi direction. 

The wave function for the vacuum state is obtained as 


lim {U\P{T) Nt \^) = (i7|0)(0|tf) ( 2 . 11 ) 

JVt->oo 


for an arbitrary gauge invariant state |T) which satisfies (0|\1/) 7 ^ 0, where P is a projection 
to the physical (gauge invariant) Hilbert space as 


P := 


n dg x Px (fix 


x£S ' 


( 2 . 12 ) 


Here E x (g x ) generates the gauge transformation at x by g x . Note that ( P ) 2 = P and 
[T, P] = 0. While we explicitly write P in the above expression since (U\ is not gauge 
invariant, the formula without P is equaly correct since PIT) = |T). 

Thus we can write 


(U\P(T) Nt \U') = 


exp 


2 Sd(U) 


rU[N T ]=U 
IU[0\=U’ 


T>U[t] exp [5c(H[t])] exp 


- 2 s d (u’) 


(2.13) 


where 


Nj 1 — 1 N t 

VU[t] : = n n dU e [t] nn dg x [t\ 

t =1 tec t=0xes 

Nt— 1 Nt 

S G (U[t ]) := E S d {U[t])+'ES 0 mt-l],U«[t]) 

t =1 t =1 


(2.14) 

(2.15) 


where U x ^[t] = gx[t]U x ^[t\g\ +jl [t\. Defining a new gauge field as U Zj 0 := gl[t\g x [t + l\ where 
z = (x, t) is a d + 1 dimensional lattice point, and introducing a new notation for gauge 
fields U Zttl with fi = 0 , 1 , ■ • • d, we have 


(U d \P(f) NT \U' d ) = 


exp 


2 s d m 


" u ZT =u d 


•Uzo =u' d 


VU / Vg ZQ exp [S p i aq (U d )} exp 


2 S d (U') 


(2.16) 
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where 


_/V T — \d d 

S P i aq (U) = iE[E E p A*) + E {-fofc( 2 o) + Pko{ z o)} 

* xGS t= 1 /-(= 0,u=0 k= 1 


(2.17) 


with zo = (x, 0) and zt = (x, IVr). Here U ZQ = Ud means U ZOt k = ( U d ) x ,k for A: = 1,2 ,--- , d. 
Note that since S'piaq and Sd do not depend on g ZQ , the gauge transformation left after the 
change of variables, we have f T>g zo = 1 in the above expression. 


Path integral expression The (unnormalized) density matrix for the vacuum state, 
can be obtain as 


(U d \p\U' d ) = e^ [Sd{Ud)+Sd{u '^ lim [ VUS{U- 

N t ^oo J z o 

\[S d {U +)+S d {U _)] 
x e t z t 


Ud)S(U 4 - U' d )e s p'^) 

(2.18) 


where i£5, — Nt <t< Nt, z q = (x, 0 ± ), and = (x, ±Nt)- 

In practice, one often employs the periodic boundary condition at ±Nt in the Euclidian 
time, which correspond to the thermal density matrix at temperature T = 1/(2 Ntcl), where 
a is the lattice spacing. In this case, after interchanging t = 0 and t = ±Nt, we have 


(U d \p T \U' d ) = e hl s d(Ud)+s d (u' d )] f 4 J VUe s^{U) 

Ju *-= U ' d 

where 

Nt — Id d 

s p ^( u ) = ^ E[ E E ^w + E{ p ofc(^)+^o(^)} 

® xes t=-N T + 1 H=0,v=0 k= 1 


(2.19) 


( 2 . 20 ) 


The density matrix for the vacuum state is reproduced from p T by the T —> 0 limit. 


Reduced density matrix for replica formula We now consider two regions V and 
V = C\V, and denote U = (Uy, Uy) and U d = ( Udv,U d y ). Then the reduced density 
matrix f)y can be written as 

Pv(Udv\U' d y) := {U d y\Ppl\U' dV ) = ^ VUe^ s ^ +s ^e s ^ {U) . (2.21) 

J U z -v=U' dV 

Z T 


The replica formula for the entanglement entropy in now given as 


S(V) = lim 


1 


n-¥ 1 1 — n 


Z„ 


T\n 


log ( ^ ) , Z n := Tr(py) 


where Z n can be expressed in the path-integral as 


Z n — 


nH 


Pv(Ui; U 2 )pl(U 2 - U 3 ) ■ ■ ■ P y(U n ; t/i) 


( 2 . 22 ) 


(2.23) 
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and py(Ui,Uj ) is given in (2.21). Note that (2.23) is invariant under the local gauge 
transformation g in d + 1 dimensions with the period 2 Nt ( not 2nNx) at the boundary, 
which satisfies 


9(x,2kN T ) = 9{x,o) at x G <9V, k = 1, 2, • • • , n. (2.24) 

3 Z N gauge theories in an arbitrary dimension 
We consider the Z^ gauge theories in this section. 

3.1 Some properties of divergence-free flux-configurations 

Flux-configuration For each link £ G C, we associate a flux ki G {0,1 ,,N — 1}. We 
assume the consistency kp = —kg. Here and throughout the present paper, equalities for 
the flux are with respect to mod N. We denote by k = (fc^eC a configuration of flux over 
the whole lattice which satisfies d%(k) = 0 at v x £ S, where 

dx( k ) ■= *(*,!/)■ ( 3 -l) 

y£S 

s.t. (a ',y)£C 

is the divergence of k at x associated with region C. We denote the set of all divergent-free 
fc’s by T. 

Take an arbitrary subset V C C. For any k G F, let Ry(k ) be the configuration 
obtained by omitting all the flux outside V. We then denote by JFy the set of k' which is 
written as k' = Ry(k) for some (not necessarily unique) k G T. 

Incoming flux and decomposition of T Fix an arbitrary subset V C C. For any 
k G JF, we define 

f v (k):={d%(k)) xeav , (3.2) 

where d\ is the divergence associated with the region V, obtained by replacing C —> V 
in (3.1). Note that fy(k) is the list of incoming flux at each site on the boundary dV. 
Recalling that dV = dV for V = C\V, we have 

fy( k ) = ~fy ( k )j (3-3) 

which represents the conservation of flux at the boundary sites. 

For a give subset V C £, we say that / G {0,1,..., N — l} a * is admissible if there 
exists at least one k G T such that fy(k ) = f. Then we have a natural decomposition 

J- = (J J-(-f), (3.4) 

/ 

where the union is over all admissible /, and 

j* f) i={k€F\f v {k) = f}. (3.5) 
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It is remarkable that all J-W with admissible f are completely isomorphic to each 
other. To see this, take arbitrary f 1 and f 2 which are admissible. Choose and fix Aq, ko G F 
such that fy(ki) = f, for i = 1,2. Then we define a map 9^2 '■ F^ 1 ^ —>• F^ 21 by 
^ 1 , 2 (fe) := k — k\ + k 2 and its inverse map 992,1 : F ^ 2 ) —>• F ( f 1 - 1 by 992.1 (fe) := fc — &2 + fei- 
Since 9 ?i, 2 (fca) / ^ 1 , 2 (fe&) for fc a / fc&, k a ^ G F^ 1 ^ and a similar property for 992 , 1 , these 
maps establish a one-to-one correspondence between the elements of 1 ) and 2 ^: F ^^ 
and F^f 21 are isomorphic to each other. 

Finally let us evaluate the number of all the admissible /’s. Decompose V and V 
into connected components as k = V\ U ■ ■ ■ U V n and V = V\ U ■ ■ ■ U V m . (For example, 
see Fig. 1 in appendix A.) Correspondingly, the boundary dV is decomposed as dV = 
dV\ U • • • U dV n = dV = dV\ U • • • U dV m . Then the divergence-free condition for k implies 
that an admissible incoming flux f = ( f x )xedv satisfies 

Y fx = °, i = 1, • • ■ ,n, 

xedVi 

Y fx = 0, j = 1, • • • ,m, 

xGdVj 

with an additional condition that 

n m 

EE'-EEi < 3 - 8 ) 

*=1 xedVi j =1 xedVj 

for an arbitrary f even without satisfying the divergent-free condition. Thus the total 
number of the admissible f ’s is readily found to be N^ 9V \ ( n+m 1 ) ; where \dV\ denotes the 
number of sites in dV. See appendix A for a more rigorous discussion. 

Decomposition of k Let V Cfbea subset, and f be an admissible incoming flux. We 
define Fy as the set of k! G Fy which is written as k' = Ry(k ) for some (not necessarily 
unique) k G F^\ 

Note that an arbitrary k G F^ is written as 

k = (ky , ky ), (3-9) 

where ky = Ry(k) and ky = Ry{k). We then have ky G Fy\ and ky G Fy ^ ■ We 
remark here that Fy J y is the set of configurations on V with incoming flux to V (i.e., 
outgoing flux from V ) equal to — f. 

3.2 Zn gauge theories 

We consider the Zjy gauge theory, generated by Zat = {g° = 1, g 1 , - ■ ■ , g^ -1 }, where g is 
a generator of the Zn and satisfies g N = 1 and g _1 = gt. 

Operators and states With each link £ G £, we associate the Wdimensional Hilbert 
space Tit, whose orthonormal bra-basis is given by i{U\ with U G Zjv- The coordinate 
operator Ug and the momentum (electric) operator Ef act on this bra-state as 

i(U\U £ = t (U\n(U), t {U\Ef = e (gU\, £ (U\Ef = e(Ugi\, (3.10) 


(3.6) 

(3.7) 
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where r\{U) is the fundamental representation of the Zn group such that r\ (< 71 ( 72 ) = 
ri(gi)r\(g 2 ) for g\,g 2 6 %n- All irreducible representations are one dimensional and ex¬ 
plicitly given by Vk(g) = e' l27rk / N for k = 0,1, 2, ■ ■ ■ , N — 1. 

The basic ket-state | h)e with h G is defined as 

e(U\h)e = 5u,h, (3-11) 


and the general state can be expressed as 


N -1 

\^)i= ^2c n \g n ) e , c n G C, 

n=0 

where |< 7 ”)^ with n G {0,1,..,,IV — 1} forms a basis of |T)^. 
We now introduce the basis of the flux representation as 


(3.12) 


N—l 

\k)i = - 7 = r k (g n )\g n ) e , k G {0, 1,..., N — 1}, (3.13) 

v n =0 

which leads to 

N- 1 TV-1 

= 7 = XI = -=r*(l7) £ 

n=0 ^IV n=0 

iV-1 

= E <(C| 3 n )l = 1. (3-14) 

n=0 

Since 

iV-l 

= -= 2 rjto")!!,”- 1 ), = rn( 9 )|k>«, (3.15) 

* n=0 

so that | k)e is an eigenstate of the electric operator E ( j with an eigenvalue r k {g)- We shall 
use this electric flux representation, which is suited for studying reduced density matrices. 
The Hilbert space 7-L for the whole system is spanned by the basis states 


|fc> :=®|A*>* 

lec. 


(3.16) 


where k G T. The gauge invariant condition at x that 

G%\k) = | k), G% := E 9 [xy) (3.17) 

yes 

S.t. ( x,y)£C 

leads to 


r 1 (g) d ^ = 1^4{k) = 0 l (3.18) 

where we use a property that r k (g) = r\(g) k . Therefore the divergence-free condition for 
k corresponds to the gauge invariance of | k) at all x G S. In terms of link variables U, 
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|fc) represents rk e (Ue) = r\{Ui) kl at each link l, and the gauge invariant (divergence-free) 
condition means that {r\(U() ke }i & c forms several closed loops with identifications that 
ri(U t ) N - k * = ri(uj) ke . 

For a subset V C L we also define My as the space spanned by 

\ky)v := (^) \k()e, (3.19) 

lav 

where ky = (k^e^y G Fy. The state \ky)y is not necessarily gauge invariant. 

»» / j ?\ ^ / j;/ \ 

If ky G Fy and k' v G Fy J with f ^ f', the corresponding kets | ky)y and | ky)y 
are orthogonal. This means that the Hilbert space My is decomposed into a direct sum 

'Hv = Q)M ( J\ (3.20) 

/ 

where My ; is spanned by \ky)y with ky G Fy . 

Fix an arbitrary subset V C £ and let V = C\V. Corresponding to the decomposition 
(3.9) of k G F^\ the state (3.16) is decomposed as 

I k) = | ky)y tg) | ky)y, (3.21) 

where \ky)y G m[P and | ky)y G My ^. 

Reduced density matrix Take an arbitrary normalized state T) G M. and expand it 
as 

|*) = X>( fc )l fc >> (3-22) 

kaF 

where ip(k) G C. We shall hx a subset V C £ and its complement V = £\V, and study 
the reduced density matrix in the region V for the state IT). 

By taking into account the decomposition (A. 10) of F, and the decompositions (3.9), 
(3.21) of k and the corresponding ket, the state (3.22) can be written as 

i*> = E E 

/ feeA(J’) 

= ^ Ipfevi ky)\ky)y <g> \ky)y, (3.23) 

f k v aF v f) k v eFf f) 

where the first sum is over admissible f. Then the corresponding density matrix is written 
as 

l*X*l=E E E ^k V ,ky)^k'y,k'y) \k V ){k'y\®\ky)(k'y\. (3.24) 

f , k v aF v s) k v eFf f) 
f k' v aP^"> k' v aPp f) 
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Since | ky) with ky 6 J-y are orthonormal, the desired reduced density matrix is readily 
found to be 


/Sy :=!*«„ [|*)(*|] 

-E E £ 'l/j(ky,ky)‘lp(ky,ky) \ky)(ky\ 
f k v &F ( v f) kyePf f) 
k' v ^y f) 

= J2pfp { v f) - ( 3 - 25 ) 

/ 

( f) 

We have here defined the density matrix on 77 (see (3.20)) by 

Pv } = — ^2 ^( k v-,ky)il){ky,ky) \k v )(ky\, (3.26) 

Pf k V £T { v f) kydF ( y f) 

k'y£p ( J ] 

where pf is obtained from the normalization condition. 

As is well-known the final expression in (3.25) implies 

S[p v ] = H[p] + ^Pf S[py ] }, (3.27) 

/ 

where H[p] = — Pf^ogpf is the (classical) Shannon entropy for the probability distri¬ 
bution of the incoming flux through the boundary. Note that the “quantum part” S[py] 
is in general obtained by diagonalizing the expression (3.26); this calculation may be non¬ 
trivial. 

It may be suggestive to observe that, in the expression (3.27), the von Neumann entropy 
S[py ] seems to reflect “intrinsic entanglement” between V and V while the Shannon 
entropy H\p\ may simply reflect the behavior of Wilson loops that touch both V and V. 

3.3 Some properties 

In the Zjy gauge theories, the density matrix can be expressed in the flux representation 
as 

P = '^,\k)Pk,k'(k' I) \k),\k') eT-Ltot, (3.28) 

k,k' 

in general 1 , where 77 tot is the full Hilbert space without gauge invariance, and tr to t p = 1 
implies J2 k Pkk = 1- The gauge invariance under the gauge transformation G 9 X and Gy 
with v x, v y implies that p kk > can be different from zero if and only if d^(k) = dy(k') = 0 
for v x, v y. This means that k and k' are divergence-free. Therefore 

tr tot p = (P\ k )Pkk'( k '\p) = ^2 ppp = tr ph P, (3.29) 

p k,k’ef P eT 


1 This form of the density matrix is more general than (3.24) for the pure state 'h). 





where tr p h represents the trace over the physical space H. 
Furthermore, the reduced density matrix is written as 


Pv-£ 22 *22 Pfov, ky', ky, k v )\kv){ky\. (3.30) 

f k v eF(f) kyCF(-f) 
k' v ef(f) 

Therefore, for \p)v G H tot y, we have 


Pv\p)v = Q (3.31) 

unless py £ Ty, so that 

tr to t,vPv = try Pv = 1, (3.32) 

where try is a trace over Hy in (3.20). In addition, we have 

f{ky\pv\ky)fi=8fji 22 p(ky,ky-,ky,ky). (3.33) 

kyef(-f') 

Therefore we can extend py in the full Hilbert space on V, Htoty without any modifications. 

The above argument shows that p and py can be regarded as the full and reduced den¬ 
sity matrices in the full Hilbert spaces without gauge constraint. The standard method then 
can be applied to prove properties of py such as positivity and strong sub-additativity[30]. 


3.4 Entanglement entropy for special states 

Factorized states and the topological state Consider a special state in which the 
coefficients in (3.22) and (3.23) factorize as 


Ipiyky , ky) — 'lpy(ky) Ipy(ky) (3.34) 

for any k = (ky, ky) £ T . Then the three summations in (3.26) can be treated indepen¬ 
dently to give 

( 22 ^v(kv)\kv)\( 22 ^v(kv)(kv |Y (3.35) 

k v £p(,P k v £p[/ ) 

which shows that py ; is pure, and hence S[py ] = 0. In this case we find that the 
entanglement entropy 5[/5y] is equal to the Shannon entropy H[p\ for the probability 
distribution of the incoming flux f. 

The topological state, in which all the coefficients ip(k) in (3.22) are identical, is an 
example where the factorization condition (3.34) is satisfied. (See Refs. [9, 10, 23, 24, 31, 32] 
for related issues.) This state is called the topological state, since an arbitrary (Wilson) 
loop has an unit eigenvalue. Namely, for £ JF, we have 

17 fc '|topo) = |topo), U k> := U(U^, (3.36) 

ter 


Pv ] - ( I'MMI 2 ) 

V -n ' 
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where 


|topo) =lf} X |fc) 
k&T 


with ^ is a complex number. Indeed, since 

l{U\U t \k) t = =r k (U)n(U) = J=r k+1 (U) = t (U\k + 1 )/, 

where we use r kl (U)r k 2 (U) = r kl+ k 2 {U), we have 

U k |topo) = if} |k + k') = i\) | k") = |topo), 

keT k"ef 


where k" = k + k' e F . 

Writing a = l^l 2 , the expression (3.35) becomes 

a 
Pf 


nopo,v 


£ 


1 


fcf/GJ, 


•(-/) 


X ( X ( ky \\ 




if) 


kyGJT, 


if) 


where 



(3.37) 


(3.38) 


(3.39) 


(3.40) 


(3.41) 


We shall argue that pf is independent of /, and hence is equal to 1/iV^' \-(n+m-i) ^ w } iere 
jy|0i |—(n+m—l) j g a num ber of independent /’s as shown in appendix A. We thus get the 
desired result 

5 , [pto P o,y] = {\dV\ - n d ) log N, n d := n + m - 1, (3.42) 

where |c?V| is a total number of boundary points, n and m are a number of disconnected 
components of V and V. respectively. 

The asserted independence is easily seen if one recalls the one-to-one correspondences 
between F^ 1 with different f. Take admissible f 1 and f 2 ■ By restricting the map ip 12 to 
Fy and Fy 1 , respectively, we obtain one-to-one correspondences between Fy 1 and 
Fy 2 * and between Fy ^ and Fy ^ 2 \ We thus find that the expression (3.40) for different 
f are in perfect one-to-one correspondences, so that a numbers of elements in both Fy ; 

~ /_ -P \ 

and Fy J does not depend on /, and thus pf is independent of f from (3.41). 

States with products of two loops We consider a simply entangled state, given by 


1 

(U |*> = -j= ( Ur v ) ki (3.43) 

vn i =1 

for n < N, where integers hi s satisfy ki 7 ^ kj for i ^ j, and Ty and Tp are closed loops 
in V and V without touching the boundary, and Ur is a product of r\(U) along the closed 
loop T. The reduced density matrix then becomes 


Pv(U,U’) 


1 n 


(3.44) 
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so that the entanglement entropy is given by 


S(V) = logn. 

In terms of the decomposition in eq.(3.27), we have 


Pf = 


1 for f = 0 
0 otherwise 


-'(/) r 'E' 

Pv ~ $f,oPvi 


so that 


H]p] = 0, PfS\p$ ] ] = 5/,o logn. 

A simply disentangled state, on the other hand, is constructed as 

<£W = l{Ur v +Up v }®{u£\+u£i }, 

which leads to 

pv(u, u') = \ [u k ' v + u k i] {( U^ + (I/f,) fc2 } , 
S(V) = 0. 


(3.45) 

(3.46) 

(3.47) 

(3.48) 

(3.49) 

(3.50) 


Single-loop states An entangled loop state is constructed as 

1 n 

(*W = - 1 =Y,^vr v ) ki (3.51) 

V i 
2—1 

for n < N, where integers k^s satisfy ki / kj for i ^ j, and IVTy is a closed loop with 
TyTy G /o / 0. The reduced density matrix and entanglement entropy are given by 

1 n 

Pv(U, U') = -J2 Up (U^ v ) k S (3.52) 

n i =l 

S( V) = logn. (3.53) 


In terms of the decomposition in eq.(3.27), we have 

_ 3 


Pf = 


1/n for f = d kif o 
0 otherwise 


Pv 


(k/o) — jjki 




\ki 


so that 


H\p\ = logn, 5[pJ/ } ] = 0. 

An example of a disentangled loop state is constructed as 

OT®> = ^ {(Cr r r t ,)*‘ + (Crvr,.)' 11 } , 

which leads to 

Pv(U,U') = \ {(U Fv ) k i + (t/ n ,) fcl } {(t/jfy)* 1 + (E^)* 1 } , 

S(V) = 0. 


(3.54) 

(3.55) 

(3.56) 

(3.57) 

(3.58) 
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3.5 One dimensional lattice without boundary 


Since one dimension is a little special, we here consider the d = 1 case separately. 

We consider Z ^-gauge theory on one dimensional lattice with periodic boundary con¬ 
dition. Note that the open boundary is incompatible with the gauge invariance. Since 
there are no Wilson loops (except one big loop on a whole lattice), only the momentum 
operator Ef is a gauge invariant operator. 

Considering the Gauss law, every link has the same electric eigenvalue. Therefore, 
physical state is given by 

l^fc) = 0 \k)i = \k) v ® \k)v (3.59) 

i 

with 

|fc)v = 0|fc),, \k)y = (g)\k) l , (3.60) 

l&v l£V 

for arbitrary partitioning. 

Topological state A topological state is given by 

! N-i 

|topo) = —= l^fc)- (3-61) 

k =o 

The global density matrix becomes 


p = |topo)(topo| 
k,k' 

= \k)v(k'\y 0 \k)y{k'\y (3.62) 

V k,k> 

and the reduced density matrix 

p v =Tr v p= ^Y\ k )v( k \ v - ( 3 - 63 ) 

k 

Therefore, the entanglement entropy of the one dimensional topological state is given by 

5(to P o)(y) = _ tT p v i og p v 

= log AT (3.64) 

= [n B - no) log N, (3.65) 


where ns is the number of boundary points and no = n + m — 1. Since n = m and n B = 2n, 
we always have 

n B — no = 1 (3.66) 

in one dimensional space. The entanglement entropy does not depend on the number of 
links in V. The result in (3.65) is the same as the topological state entropy formula in 
d > 2 lattice, 
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General state We consider general state as 


N—l 


\a) = 


Kj2 a ( k )\*k) 


k =0 


with the normalization coefficient 


K = 


Ek=oH k )\ 2 ' 

The reduced density matrix is given by 

pv = K 2 ^ a(k)a(k)\k)v(k\v■ 

k 

The entanglement entropy is given by 

5(1/) = -trpy \ogp v = - ^p fc logpfc 
with pk = \oi(k)\ 2 /Yh \ a ( k )\ 2 - For the topological state po = pi = • • • = Pn-i = 1/-/V, 

5(to P o)(y) = iogj\r. 

For pure state po = l,pi = • • ■ = Pn- i = 0, 

5(pure) (y) = 0 . 

A simply entangled state 


with k A k !, gives 


l*> — {l^)v ® | k )v + \ k ')v <8> I k ')v} 


S(V) = log 2. 


(3.67) 

(3.68) 

(3.69) 

(3.70) 

> 

(3.71) 

(3.72) 

(3.73) 

(3.74) 


3.6 Relation to other proposals 

We here discuss relations of our definition of entanglement entropy (or the reduce density 
matrix) for gauge theories, in particular, the Z jv gauge theory to other proposals. 

Our definition is equivalent to the electric boundary condition (electric center) in Ref. [23] 
and in Ref. [24], to the extension of the Hilbert space in Ref. [25], and to the extended 
lattice construction in Ref. [26]. In this definition, the reduce density matrix py, from the 
whole density matrix p restricted to the region V, satisfies 

( Oy ) := tr [Oyp\ = tr y[O v pv} (3.75) 

for v 0y E Ay, where Ay is the set of gauge invariant operators on V, generated by 
with l E V and U p with the plaquette whose links are all included in V, and try is the 
trace over kLy. It is noted that Ay is the maximal gauge invariant algebra on V. 
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The trivial center definition in Ref. [23], denoted by p ( y, is equivalent to the gauge fixed 
theory where the boundary links in the maximal tree are all fixed to the unit element. 
In this case, however, the set of gauge invariant operators Ay, generated by Ef with 
£ 6 U\{maximal tree} and the same set of plaquette U p on V, is smaller than Ay. Similarly, 
the algebra Ay associated with the magnetic center[23, 24] is smaller than Ay . Therefore 
both Ay and Ay do not represent the region V algebraically, so that definitions based on 
the trivial center and the magnetic center are inadequate for the entanglement entropy or 
the reduced density matrix on the region V. 

In conclusion, our definition of the entanglement entropy or reduced density matrix 
gives the unique definition of these quantities on the region V, in the sense that our reduced 
density matrix is associated with the maximally gauge invariant algebra Ay on V. 

3.7 Gauge fixing 

Since the reduced density matrix py does not have the full gauge invariance as mentioned 
before, the entanglement entropy may depend on whether gauge fixing is employed or 
not in the calculation, and on the choice of the gauge if the gauge fixing is used. In 
this subsection, using a simple example, we explicitly demonstrate that the entanglement 
entropy with some gauge fixing is different from the one calculated without gauge fixing. 

We consider the Zjy gauge theories in one dimension with periodic boundary condition 
in subsection 3.5. Without gauge fixing, the entanglement entropy is given in (3.70) as 

S{V) = -VVfogPfc, Pk = K 2 \a(k)\ 2 , K 2 = } (u \\2 ( 3 -76) 

for a general state 

\a)=Kj2<*(k)\*k)- (3-77) 

k 

Take L lattice points on the circle as S = {1, 2, • • • , L} and C = {(1, 2), (2, 3), • • • , (L, 1)}. 
Links in the region V are given by Cy = {(1, 2), (2, 3), • • • , {Ly — 1, Ly)}, while those in V 
by Cy = {{Ly, Ly + 1), • • • , (L, 1)}, where 0 < Ly < L and dV = {1, Ly}. Using gauge 
transformations on all points in S except one, we can always make U( = 1 for all £ £ C 
except one £ which may be in Cy or Cy. In any cases, the reduced density matrix from 
the global pure state is always pure, so that the entanglement entropy is always zero. This 
is clearly different from (3.76) without gauge fixing. 

We next consider the gauge fixing using all points in S except dV. In this case we can 
make U(_ = 1 for all t E C except two Us, one £ in Cy and the other in Cy. For example, 
we can take U( 12 ) A 1 and LW i) A 1- Since the gauge invariance still holds on the site 1, 
the physical state can be written as 

I a) = ^5>(A0l%- 2 ) ® \ k ){L,i). (3.78) 

k 

Then the reduce density matrix is given by 

py = K 2 Y,Hk)\ 2 \k) (12) (12) (fc|, (3.79) 

k 
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which leads to (3.76) for the entanglement entropy. For the topological state, it reducers 
to 

S {topo \V) = \ogN. (3.80) 

The above consideration leads to an important lesson that the entanglement entropy 
does not depend on the gauge fixing if and only if points in dV are excluded in the gauge 
fixing (including no gauge fixing at all). Otherwise, the entanglement entropy does depend 
on the gauge choice. 

4 Conclusion 

We have proposed the definition of the entanglement entropy in lattice gauge theories for 
an arbitrary subset of links not only in abelian theories but also in non-abelian theories, 
and explicitly given the replica formula based on our definition. In the Zjy gauge theories, 
we have expressed the whole Hilbert space by the flux representation basis states which are 
eigenstates of the gauge transformations. By using these basis states, we have explicitly 
argued that all the standard properties of entanglement entropy hold in our definition and 
calculated the entanglement entropy for topological states as 

S[p topo ,V} = (\dV\-n a )\ogN. (4.1) 

We have also found that the entanglement entropy depends on the gauge fixing in general. 

It will be important to extend our analysis for the Zn gauge theories to non-abelian 
gauge theories, since our definition is applicable also to non-abelian cases without any 
difficulties. In order to calculate the entanglement entropy analytically in non-abelian gauge 
theories, we need some useful basis such as the flux representation in the Zn gauge theories. 
In the Zjsr gauge theories, the flux representation basis diagonalizes gauge transformations 
simultaneously. On the other hand, in non-abelian gauge theories, gauge transformations 
cannot be diagonalized simultaneously since they do not commute each other. We therefore 
need some new ideas for non-abelian gauge theories. In appendix B, some analyses in this 
direction are given. For example, the entanglement entropy for the topological state in one 
dimension is calculated as 

Sv = log |G| (4.2) 

in the discrete non-abelian gauge theories, where |G| is a number of elements of the discrete 
group. 

Others directions in future investigations include perturbative calculations for the en¬ 
tanglement entropy in gauge theories [33-36] without gauge fixing at boundaries and nu¬ 
merical simulations for the entanglement entropy in lattice gauge theories[25, 38, 39]. 


After completing our investigations presented in this report, we noticed a paper [40] in 
which the authors also propose the definition of the entanglement entropy in lattice gauge 
theories. We find that their proposal is identical to ours, though research directions in this 
paper are somewhat different from theirs. See also Ref. [41] for a related result. 
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A The number of admissible f 

Here the whole L is assumed to be finite and connected. 

Suppose that V and V are decomposed into connected components as 

V = V 1 U---UV n , V = V 1 U---UV m , (A.l) 

with n, m > 1. Consequently the boundary dV = dV is decomposed as 

8V = dVi U • • • U dV n , dV = dV x U • • • U dV m , (A.2) 

where dVi and dVj are the boundaries of Vi and Vj, respectively; they may not be connected. 
Let us denote by 

?-={f = (fx)xedV I f x € {0,1,..., N - 1} for s € dV } (A.3) 

the set of all configurations of incoming currents (including “unphysical” ones). We have 

\T\ = xi av 'L 

The necessary and sufficient conditions for the admissibility of f are 

Y fx = 0 , and Y I* = °> ( A - 4 ) 

xedVi xedVj 

for all i = 1 ,n and j = 1,m. There are n + m constraints, but they are not 
independent. To see this note that any f satisfies 

n m 

EEt=EEi' < A ' 5 ) 

i =1 xGdVi j =1 x&dVj 

because dV = dV. We thus see that 

n m —1 

E/- = EEt-EEt < A - 6 ) 

xedVm i=1 xedVi j =l xedVj 

holds for any f. It is therefore sufficient consider the constraints (A.4) for i = 1,..., n and 
j = 1 ,,m — 1. There are 7 := n + m — 1 constraints. 
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To count a number of admissible /, let us introduce matter fields (or external sources) 
with Zn charge on lattice sites {rrjxes- For a given charge density distribution {g(x)} xg5 , 
an admissible flux in this general case is determined so as to satisfy the Gauss law as 


£ fx = QiS{0,l,...,N-l}, (A.7) 

x£dVi 

for each i = 1 ,..., n, and 

J2 f* = -Qo+n € {0,1,..., N - 1}, (A.8 ) 

x&dVj 

for each j = 1 ,... ,m, where Qk (k = 1, • • • ,7 + 1), a sum of q(x) over inner point, is a 
total charge inside the region Vi or Vj excluding boundaries. The minus sign in the second 
equation comes from the fact that a flux on dV has a relative minus sign with respect to 
a flux on dV. Due to the constraint (A. 6), we have 

7 

Q'y+i = - y ^Qk, (a.9) 

k=l 

so that only Q 1, • • • ,Qy are independent. We then define as the set of f E T 

which satisfies (A.7) and (A. 8). It is then easy to see 

F = U •Fqi,...,<3 7 - (A.10) 


Note that To.-.-.o is the set of admissible /’s that we are interested in. 

Now we will argue that J- o,... ,o is isomorphic to J r Q lt ... n for an arbitrary Q 1, • • • , Qy. 
Take one internal point xj~ from each region V) or Vj. Connect these points by the following 
condition: (1) links can be used once. (2) except start and end points, each point belongs 
to only two links (3) the end point is always It is easy to see such a connection 

always exist. By changing the order of point along this connection and renaming xk in 
this order, we write the connection as T XlX2 T X2X3 ■ ■ ■ T Xi _ lX ^T XlX ^ +1 , where T x% xi+\ is a set 
of links which connect Xi and Xi+ \. For an illustration, see Fig. 1. 

For v {Qi, • • • , Q 7 } (this is also reordered), we define k® 1 ’"' ’^ 7 on a link i as 


k 


Ql,-" ,Q~t _ | Si=l Qii f° r £ G Fx k x k+1 
I 0, otherwise 


(A.ll) 


See Fig. 1 again as an example. A blue letter such as Q 1 + Q 2 represents a charge on some 
lines, while a red letter such as Qt~ is a charge on the point x Note that the net charge 
flowing out from the k -th region (some Vi or Vj ) is equal to Yli= 1 Qi ~ Yli=i Qi = Qk- It 
is then easy to see that the map for k E To,--- ,0 defined by 

^Po- ’-'o' Ql ( k ) -=k + k Ql ’-’ Q i (A. 12) 


establishes an isomorphism from To,-- ,0 to Tq^... t Q . This proves the number of t Q 

is independent of Qi, ■ ■ ■ , Q 7 . 

A number of possible charge distribution {Q\, Q 2, • • • , Q-y} is N 7 = j\r n+m_1 . There¬ 
fore, for any charge distributions {(/(aOjies including {(7(x)}xe*s = {0} x e*S) the total num¬ 
ber of the admissible f is ]\f\9V\-(n+m-i) _ 
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Figure 1. An example of the connection and charge distributions. 

B Entanglement entropy for non-abelian gauge theories 
B.l About the Hilbert space on a link 

We generalize the formulation of the Z n case to non-abelian gauge theories. We take a 
group G which we assume to be a compact group. We define the momentum operator 
Lf(g) and the position operator UJ via 

{U\L e (g)\V) = Vig-'U), {U\(Uj) a ^) = tt (U) a ^(U) (B.l) 

where g G G and tt is a representation of G. If we inverse the direction of link £. the 
operator L i r(g) and UJ T is defined as follows: 

(U\L lT (g)\*) = *(Ug), (U\(U£ T ) a p\V) = tt (CT 1 ) 0 ^ (17) = (tt(U)*)^(U) (B.2) 

It is known that the L 2 space on a group G (square integrable functions over G) 
decomposes to the direct sum of tt^ IE tt which is a irreducible representation of G x G as 
follows [42]: 

L \G)^ © Vj ® 14, (B.3) 

Trelrr(G) 

where we denote tt as an (unitary) irreducible representation of G and Irr(G') as the set of 
irreducible representation and tt' (g) = t n(g^ 1 ) is the dual representation. The meaning of 
(B.3) will become clear below. 

We first consider the basic state \TT a g) defined via 

(U\TT a g) = TT\U)g a , (B.4) 

with which we can explicitly write the action of Li(g) as 

(UlLiig)^) = tT\ g- l U)f = ^{U^gTg = ir{g)'^(U\ tt%). (B.5) 
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Therefore we have 


Le{g)\Tr a p) = \'K a 1 )i{(g)' 1 p, (B.6) 

L er(g)\v a p) = n{g- l )%Wp). (B.7) 

The dual vector of |7r“ ( g)is given by (|7r“^))^ = (7r^ Q |, and the projection operator to 
the subspace 7U Kl n is given by 

Pn = diml4^ |7T Q /3 )(7r /3 Q |. (B.8) 

a,(3 

The factor dim V n is needed here since the normalization of vector y a p) is given by 

KpWs) = J du ^piUWsiU- 1 ) = (B.9) 

By the projection (B.8), the meaning of (B.3) becomes clear. 


B.2 Gauge invariant states 

In the lattice gauge theory, the total Hilbert space Hq is ^>i(L 2 (G))i. The physical Hilbert 
space T~i as the subspace of T~Lq is consist of gauge invariant states, which satisfy 

om = I*), G° x := n L t (g), (B.10) 

yes 

s.t. e=(x,y)ec 

at v x E S. The basis state in Tio is written in general as 

I *“„>:= 0lfa)“‘ft> (B. 11) 

tec 

where n£ indicates an irreducible representation of G on a link Z. 

Unlike the Zn gauge theories, it is not so easy to write gauge invariant conditions for 
the state in (B.ll). Let us consider the one dimensional case as a simplest example. In 
this case, the nontrivial part of the gauge invariant condition at x becomes 


7r ai 7i)u 7r (5') 7l /3i 1 )“ 2 72 W l2 f3 2 )t 2 = k ai pi)ii ® K“ 2 / 3 2 )i 2 (B.12) 


where l\ = {x — l,x) = (x,x — 1) T and Z 2 = {x,x + 1). Integrating this equation over g 
with f dg = 1, we find that a gauge invariant state at x has a form as 


1 

dim 14 


7r a 7 )h ® |t Gp)e 2 , 


(B.13) 


where two irreducible representations on Z\ and Z 2 must be equal. 

In higher dimensions, however, the condition becomes more complicated. On a d- 
dimensional hyper-cubic lattice, the gauge invariant condition at x reads 





7 r 


- 'IPr 


(B.14) 
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where i M = (x,x + n) and £p, = (x,x — fi). This implies that a product of 2d irreducible 
representations of tt IJ: and 7 must contain the trivial representation. For example, in the 
case of SU(2) gauge group at d = 2, 4 non-negative integers £ 1 , 2 , 3 , 4 , which are numbers of 
boxes in the SU(2) Young tableaux and specify irreducible representations of SU(2), must 
satisfy 


{|£i - k 2 \, ■■■ ,£1 + £2} n {|£ 3 - fc 4 1 , • • • ,£3 + £4} / 0 - 

For general gauge groups in higher dimension, it is hard to find a simple condition for 
(B.14). 

B.3 Examples 

As was seen in the previous subsection, it is not so easy to construct general gauge invariant 
states in higher dimensions. Therefore, in this subsection, we consider two examples at 
d= 1. 

B.3.1 One dimensional topological state with periodic boundary condition 

Assume that there are N links on a circle(i.e. the periodic boundary condition). In this 
boundary condition, similar results are obtained by Donnelly[27] in the theories defined on 
the continuum space. The physical Hilbert space is given by the gauge invariant functions. 
From the analysis in the previous subsection, the basis are given by the characters of 
irreducible representations as 

H>)= ^2 V’OOK), | if)) eu 

7r£lrr (G) 

|7r) = |7r Ql a 2 ) ® K" 2 a 3 ) ® ‘ ‘ ‘ ® K^ai), (^K) = 1- 
The value of 17r) at \U\, • • • , C/jv) = \Ui) 0 • • • (8) |L£v) becomes as follows. 

(Ui,--- ,U N \n) =7T\U 1 ) a2 a ^\U 2 ) a3 a2 ---n\U N ) ai aN 

= tr(7r t (t/ JV - ■■[/!))• (B.17) 

As we have done in the abelian cases, to divide the physical Hilbert space into the tensor 
product of Hilbert spaces on the region V and V, we embed the physical Hilbert space T~i 
into a larger Hilbert space T~L' where 2 

«' = © © (B.18) 

7r,7r / ma,m(3,rna'rri(3' 

Here ma = (ai, • • • , gm), m/3 = (j3\, ■ • • , /3m) are labels of boundaries when the subsystem 
V is consist of M intervals and ma ', m(3' are the corresponding ones in V. Then we trace 
over regarding the physical wave function |V>) as an element of . 

2 Unlike the Zn gauge theories in the main text, we here consider the minimum extension where gauge 
invariance is abandoned only at boundaries. 


(B.15) 

(B.16) 
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As the simplest case, we consider V (and V) is an interval. In this case, the basic is 
written as 

|tt) = (dimK-) -1 x ^ VdimVnlirvp) <8> VdimI4|vr^ Q ). (B.19) 

Oiy(3 

The reduced density matrix for physical wave function (B.15) is given by 

pv= p(7r)(dim 14-) 2 (■\J dim K-1 -Kyp) )(y 7 dim V n (TT ^ a [) (B.20) 

-n-elrr(G) a,/3 

where p{? r) = | , 0(7r)| 2 . Its entanglement entropy is given by 

S V = ~ ^2 P( 7r )(dimI4)“ 2 log(p(7r)(dimI4)“ 2 ) 

7r£lrr(G),a,/3 

= ^2 P( 7r )( — l°gP( 7r ) + 2 log dim 14) (B.21) 

TrGlrr(G) 

Using the above result, we compute an entanglement entropy of the topological state 
in finite non-abelian group G. The topological state is given by 

^ _ 

|t°po> = /\n\-m+i £ ^ 9 n) ® Ig^gi) ® • ® Ig^dN-i), (B.22) 

vl G r 9i£G 

where |G| is the number of the element of G, and states satisfy (g\h) = \G\5 g> h- Here |topo) 
is written as the element of Ho, though it is gauge invariant, and the coefficient ip(n) is 
given by 

V’M = (vrltopo) = (B.23) 

which leads to p(n) = (dim I4) 2 /|G|. Thus the entanglement entropy is calculated as 

Sv= E (^W!(-k«a^£ + 21o*dtaV.)-k« l G|, (B.24) 

where we use the identity X^7reirr(G)(di m I4) 2 = |G|. This result agrees with (3.64) for the 
Z]y gauge theories 

B.3.2 One dimensional topological state with open boundary condition 

Next we consider the case with open boundary condition. From the gauge invariance of 
the bulk, physical wave functions are given by a linear combination of functions on G as 

W) = ^2 V’( 7r )«^ V dim 141 ntot a p) IV’) £ % (B.25) 

7rGlrr(G),Q,/3 

|vr to t a /3) = ® k Q1 Q2 ) (8) ■ ■ ■ (8) \^ aN ~ 1 g) (B.26) 

The value at \U\, ■ ■ ■ , Un) becomes as follows. 

(Ui, • • • , U N \7T to t a g) = 7r t (Hi) ai “7rt(t/ 2 ) Q2 ai • • • tt\U n )^ 

= ^{U N ---U 1 ) p a . (B.27) 
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This confirms that the physical Hilbert space is spanned by the functions on the group G. 

For example, we consider V is an interval in the middle. In this case, the basis is given 
by 

VdimK-lvrtot 0 /?} = (dimT4) _1 ^ i/dim V^|7T^ 7 ) <g> v / dimT4|7r^ 5 ) <g> ^/dirn 

7,(5 

(B.28) 

From the decomposition, we find the reduced density matrix is given by 

Pv = y^p(7r)(dimI4)~ 2 ^(y^im K-|7r^))(\/dim K-(7rf 7 |) (B.29) 

7T 7,(5 

where p{i r) = ^ V ; ( 7r )a^V’*( 7r )/3 a - The expression of the reduced density matrix is the same 
with the case of periodic boundary condition (B.20), so that the entanglement entropy is 
given by the same formula (B.21) . 

The topological state with open boundary condition is given by 

l top °) = /, -i ® l&V) ® ■ ■ ■ ® \9n-i9n-2) ® bw-i)- (B.30) 

Vl G l g^G 

Thus VdimI4-|7rtot a p) component is obtained as 


V ; a /3 ( 7r ) = v / dimT4(vr /3 a |topo) 



(B.31) 


and p(n) becomes 


p(n) 


(dim 14) 2 

\G\ 


(B.32) 


which is identical to the result with the periodic boundary condition case. We thus obtain 
the same result also for the entanglement entropy as 


S V = log |G|- 


(B.33) 


References 

[ 1 ] S. Ryu and T. Takayanagi, “Aspects of Holographic Entanglement Entropy,” JHEP 0608 , 
045 (2006) [hep-th/0605073]. 

[2] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from 
AdS/CFT,” Phys. Rev. Lett. 96 , 181602 (2006) [hep-th/0603001]; V. E. Hubeny, 

M. Rangamani and T. Takayanagi, “A Covariant holographic entanglement entropy 
proposal,” JHEP 0707 (2007) 062 [arXiv:0705.0016 [hep-th]]; T. Nishioka, S. Ryu and 

T. Takayanagi, “Holographic Entanglement Entropy: An Overview,” J. Phys. A 42 (2009) 
504008; T. Takayanagi, “Entanglement Entropy from a Holographic Viewpoint,” Class. 
Quant. Grav. 29 (2012) 153001 [arXiv: 1204.2450 [gr-qc]]. 

[3] B. Swingle, “Entanglement Renormalization and Holography,” Phys. Rev. D 86, 065007 
(2012) [arXiv:0905.1317 [cond-mat.str-el]]. 


- 25 - 















[4] B. Swingle, “Constructing holographic spacetimes using entanglement renormalization,” 
arXiv:hep-th/1209.3304 

[5] M. Van Raamsdonk, “Building up spacetime with quantum entanglement,” Gen. Rel. Grav. 
42, 2323 (2010) [Int. J. Mod. Phys. D 19, 2429 (2010)] [arXiv: 1005.3035 [hep-th]]; M. Van 
Raamsdonk, “Comments on quantum gravity and entanglement,” arXiv:0907.2939 [hep-th]. 

[6] M. Nozaki, S. Ryu and T. Takayanagi, “Holographic Geometry of Entanglement 
Renormalization in Quantum Field Theories,” JHEP 1210, 193 (2012) [arXiv:1208.3469 
[hep-th]]; A. Mollabashi, M. Nozaki, S. Ryu and T. Takayanagi, “Holographic Geometry of 
cMERA for Quantum Quenches and Finite Temperature,” JHEP 1403, 098 (2014) 

[arXiv: 1311.6095 [hep-th]]; M. Miyaji, S. Ryu, T. Takayanagi and X. Wen, “Boundary States 
as Holographic Duals of Trivial Spacetimes,” arXiv:hep-th/1412.6226 

[7] M. Nozaki, T. Numasawa, A. Prudenziati and T. Takayanagi, “Dynamics of Entanglement 
Entropy from Einstein Equation,” Phys. Rev. D 88, 026012 (2013) [arXiv:1304.7100 
[hep-th]]; J. Bhattacharya and T. Takayanagi, “Entropic Counterpart of Perturbative 
Einstein Equation,” JHEP 1310, 219 (2013) [arXiv:1308.3792 [hep-th]]. 

[8] T. Faulkner, M. Guica, T. Hartman, R. C. Myers and M. Van Raamsdonk, “Gravitation 
from Entanglement in Holographic CFTs,” JHEP 1403, 051 (2014) [arXiv:1312.7856 
[hep-th]]; N. Lashkari, M. B. McDermott and M. Van Raamsdonk, “Gravitational Dynamics 
From Entanglement Thermodynamics,” JHEP 1404, 195 (2014) [arXiv:1308.3716 [hep-th]]. 

[9] M. Levin and X. G. Wen, “Detecting Topological Order in a Ground State Wave Function,” 
Phys. Rev. Lett. 96, 110405 (2006) [arXiv:cond-mat/0510613]. 

[10] A. Kitaev and J. Preskill, “Topological entanglement entropy,” Phys. Rev. Lett. 96, 110404 
(2006) [arXiv:hep-th/0510092]. 

[11] B. Hsu, M. Mulligan, E. Fradkin and E.A. Kim, “Universal entanglement entropy in 2D 
conformal quantum critical points,” Phys. Rev. B 79, 115421 (2009) [arXiv:0812.0203]. 

[12] H. Li and F. D. M. Haldane, “Entanglement Spectrum as a Generalization of Entanglement 
Entropy: Identification of Topological Order in Non-Abelian Fractional Quantum Hall Effect 
States,” Phys. Rev. Lett. 101, 010504 (2008) [arXiv:0805.0332 [cond-mat.mes-hall]]. 

[13] S. T. Flammia, A. Hamma, T. L. Hughes, and X.-G. Wen, “Topological Entanglement Renyi 
Entropy and Reduced Density Matrix Structure,” Phys. Rev. Lett. 103, 261601 (2009) 
[arXiv:0909.3305 [cond-mat.str-el]]. 

[14] M. B. Hastings, I. Gonzalez, A. B. Kallin, and R. G. Melko, “Measuring Renyi Entanglement 
Entropy in Quantum Monte Carlo Simulations,” Phys. Rev. Lett. 104, 157201 (2010) 
[arXiv:1001.2335 [cond-mat.str-el]]. 

[15] L. Susskind and J. Uglum, “Black hole entropy in canonical quantum gravity and superstring 
theory,” Phys. Rev. D 50, 2700 (1994) [hep-th/9401070]. 

[16] D. N. Rabat, “Black hole entropy and entropy of entanglement,” Nucl. Phys. B 453, 281 
(1995) [hep-th/9503016], 

[17] D. N. Rabat and M. J. Strassler, “A Comment on entropy and area,” Phys. Lett. B 329, 46 
(1994) [hep-th/9401125], 

[18] N. Shiba, “Entanglement Entropy of Two Black Holes and Entanglement Entropic Force,” 
Phys. Rev. D 83, 065002 (2011) [arXiv:1011.3760 [hep-th]]. 


- 26 - 


[19] N. Shiba, “Entanglement Entropy of Two Spheres,” JHEP 1207, 100 (2012) 

[arXiv: 1201.4865 [hep-th]]. 

[20] T. Nishioka and T. Takayanagi, “AdS Bubbles, Entropy and Closed String Tachyons,” JHEP 
0701, 090 (2007) [hep-th/0611035], 

[21] I. R. Klebanov, D. Kutasov and A. Murugan, “Entanglement as a probe of confinement,” 
Nucl. Phys. B 796, 274 (2008) [arXiv:0709.2140 [hep-th]]. 

[22] A. Lewkowycz, “Holographic Entanglement Entropy and Confinement,” JHEP 1205, 032 
(2012) [arXiv:1204.0588 [hep-th]]. 

[23] H. Casini, M. Huerta and J. A. Rosabal, “Remarks on entanglement entropy for gauge 
fields”, arXiv:hep-th/1312.1183. 

[24] D. Radicevic, “Note on Entanglement in Abelian Gauge Theories”, arXiv:hep-th/1404.1391. 

[25] P. V. Buividovich and M. I. Polikarpov, “Entanglement entropy in gauge theories and 
holographic principle for electric strings,” Phys, Lett. B670, 141 (2008). 

[26] W. Donnelly, “Decomposition of entanglement entropy in lattice gauge theories”, Phys. Rev. 
D85, 085004 (2012). 

[27] W. Donnelly, “Entanglement entropy and nonabelian gauge symmetry”, 
arXiv:hep-th/1406.7304. 

[28] M. Creutz, “Gauge fixing, the transfer matrix and confinement on a lattice”, Phys. Rev. 

D15, 1128 (1977). 

[29] M. Lusher, “Construction of a selfadjoint, strictly positive transfer matrix for Euclidean 
lattice gauge theories,” Comm. Math. Phys. 54, 283(1977). 

[30] See, for example: M. A. Nielsen and I. L. Chuang, “Quantum computation and quantum 
information,” Cambridge University Press, Cambridge, 2000. 

[31] A. Hamma, R. Ionicioiu, P. Zanardi, “Ground state entanglement and geometric entropy in 
the Kitaev’s model,” Phys. Lett. A 337, 22 (2005). 

[32] A. Hamma, R. Ionicioiu, P. Zanardi, “Biparticle entanglement and entropic boundary law in 
lattice spin systems,” Phys. Rev. A 71, 022315 (2005). 

[33] T. Nishioka, “Relevant Perturbation of Entanglement Entropy and Stationarity,” Phys. Rev. 
D 90, no. 4, 045006 (2014) [arXiv: 1405.3650 [lrep-th]]. 

[34] V. Rosenhaus and M. Smolkin, “Entanglement Entropy: A Perturbative Calculation,” JHEP 
1412, 179 (2014) [arXiv: 1403.3733 [hep-th]]. 

[35] V. Rosenhaus and M. Smolkin, “Entanglement Entropy for Relevant and Geometric 
Perturb at ions,” arXiv: hep-th/1410.6530. 

[36] K. W. Huang, “Central Charge and Entangled Gauge Fields,” arXiv:1412.2730 [hep-th]. 

[37] P. V. Buividovich and M. I. Polikarpov, “Numerical study of entanglement entropy in SU(2) 
lattice gauge theory,” Nucl. Phys. B 802, 458 (2008) [arXiv:0802.4247 [hep-lat]]. 

[38] Y. Nakagawa, A. Nakamura, S. Motoki and V. I. Zakharov, “Entanglement entropy of SU(3) 
Yang-Mills theory”, PoS LAT 2009, 188 (2009) [arXiv:0911.2596 [hep-lat]]. 

[39] Y. Nakagawa, A. Nakamura, S. Motoki, V.I. Zakharov, “Quantum entanglement in SU(3) 
lattice Yang-Mills theory at zero and finite temperatures,” PoS Lattice2010, 281 (2010). 


- 27 - 


[40] S. Ghosh, R.M. Soni, S.P. Trivedi, “On the entanglement entropy for gauge theories,” 
arXiv: hep-th / 1501.02593 . 

[41] L. Y. Hung and Y. Wan, “Revisiting Entanglement Entropy of Lattice Gauge Theories,” 
arXiv: hep-th/1501.04389. 

[42] R. W. Carter, I. G. MacDonald and G. B. Segal, in “Lectures on Lie groups and Lie 
algebras,” Cambridge University Press, 1995 


- 28 - 


